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ABSTRACT 

Context. It is generally assumed that a large fraction of stars are initially born in clusters. However, a large fraction of 
these disrupt on short timescales and the stars end up belonging to the field. Understanding this process is of paramount 
importance if we wish to constrain the star formation histories of external galaxies using star clusters. 
Aims. We attempt to understand the relation between field stars and star clusters by simultaneously studying both in 
a number of nearby galaxies. 

Methods. As a pilot study, we present results for the late- type spiral NGC 4395 using HST/ACS and IIST/WFPC2 
images. Different detection criteria were used to distinguish point sources (star candidates) and extended objects (star 
cluster candidates). Using a synthetic CMD method, we estimated the star formation history. Using simple stellar 
population model fitting, we calculated the mass and age of the cluster candidates. 

Results. The field star formation rate appears to have been roughly constant, or to have possibly increased by up to 
about a factor of two, for ages younger than ^300 Myr within the fields covered by our data. Our data do not allow 
us to constrain the star formation histories at older ages. We identify a small number of clusters in both fields. Neither 
massive (> 10^ Mq) clusters nor clusters with ages > 1 Gyr were found in the galaxy and we found few clusters older 
than 100 Myr. 

Conclusions. Based on our direct comparison of field stars and clusters in NGC 4395, we estimate the ratio of star 
formation rate in clusters that survive for 10^ to 10** years to the total star formation to be F ~ 0.03. We suggest that 
this relatively low F value is caused by the low star formation rate of NGC 4395. 
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1. Introduction 

] It is commonly assumed that most (if not all) stars are 
■ formed in clusters. Clusters, due to dynamical and stellar 
' evolution, dissolve (|Spitzeijll987() and the stars that belong 
, to them become part of the field stellar population of the 
galaxy. To use clusters as effective tools to constrain the 
[ formation and evolution of galaxies, it is necessary to im- 
. prove our understanding of what fraction of stars end up 
' as members of bound clusters and in the field, respectively. 
Both field stars and star clusters have been studied 
extensively in the Milky Way and nearby galaxies. In 
principle, field stars hold information about star forma- 
I tion histories over the entire Hubble time and can there- 
fore provide important information abou t : gala xy forma- 
tion and evolution, e.g.. lEdvardsson et al.l |1993). Main se- 
quence stars are present at all ages, and observations of 
the main sequence turn-off can constrain epochs of stel- 
lar formation especially in dwarf galaxies where relatively 
distinct bursts are often observed (Mateo 1998). Other fea- 
tures of the color- magnitude diagram (CMD), such as the 
sub-giant and horizontal branches, giant and asymptotic 
branches, red clump stars, and red and blue supergiants 
can provide information about specific epochs of star for- 
mation. However, it is necessary to apply a more sophis- 
ticated modeling of the CMD to reconstruct star forma- 
tion histories, taking into account different effects, e.g.. 



incompleteness, resolution, depth of the observations , ex- 
tinction and chemical composition. iTosi et ahl ()1991l) de- 
veloped a method that takes these effects into account and 
attempts to use all the information available in a CMD 
to reconstruct the field star formation history of a galaxy. 
The synthetic CMD method creates a synthetic popula- 
tion that is compared to the observed CMD to constrain 
the star formation history (SFH) of the field stars in a 
galaxy. Many subsequent studies have refin ed this method 
(Dolphin 1997; Harris & Zaritskv 2001; Skill man fc GallartI 
|2002[ ). The SFHs of many galaxies in the Local Group 
and near by have been studi ed using t he synthetic CMD 
metho d: iBrown et aP (I2008D for M31 iHarris fc Zaritskv! 
(^2009*) for the LMC. IHarris fc Z aritek^ (20041 for the SMC. 
Barker j;t al. (20 07|) fo r the M33, Cole et al. ( 20071) f or Leo 
A, I Young et all (I2007D for Phoenix lAnnibaU et al.l (I2009D 
forNGC 1705. 1 Williams eTall (l2009[) f or M Sl.lLarsen et al.l 
(I2007h for NGC1313 and iReikuba etall (|2004D for NGC 
5128, among many others. 

On the other hand, the study of cluster systems 
and disruption processes can provide important insight 
into the origin of field stars in a galaxy. Clusters dis- 
rupt by means of a variety of mechanisms, includ- 
ing "infant mortality" , stellar evolution, two-body relax- 
ation, and tidal shocks that have b een extensively stud - 
ied bv iBoutloukos fc Lamerj (|2003D : iLada fc Ladal (|2003D : 
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(2008) among others. The analysis of star clusters is of- 
ten based on a comparison between the observed spectral 
energy distribution and theoretical models, which provides 
information about the ages and masses of the studied clus- 
ters and, in turn, their dissolution. However, the compar- 
ison of models with observations also remains affected by 
many uncertainties, e.g., binarity, completeness effects. 

Some studies have started to address the relation be- 
tween field stars and clusters more explicitly. The num- 
ber of stars that were formed in clusters has been es- 
timated to be 70% — 90% in the solar neighbourhood, 
while 50% — 95% of these embedded clusters dissolve i n 
a few Myrs ("Lada fc Ladal 120031: E amers fc Gielesl[2008h . 
iGieles fc Ba stian (2008) estimated that only 2-4% of the 
global star formation rate in the SMC happened in bound 
star clusters. It is currently unknown what fraction of stars 
are initially born in clusters in SMC, but if this fraction 
were as large as in the Solar neighbourhood this would im- 
ply t hat there is als o a large infant mortality rate in the 
SMC. lBastianI (|2008[ ) studied the relation between the clus- 
ter formation rate and the star formation rate (F = g^^ ) 
using archival Hubble images for high star formation rate 
galaxies and additional galaxies/clusters from literature. He 
found this fraction to be F ~ 0.08 and thus concluded that 
the fraction of stars formed in (bound) clusters rep resents 
8% of the total star formation. On the other hand, iGieled 
(|2009t) found that F is given by 0.05 < F < 0.18 in the 
galaxies M74, MlOl, and M51. In most of these cases, how- 
ever, it is difficult to tell whether there is a genuine "field" 
mode of cluster formation, or whether all stars form ini- 
tially in clusters of which a large fraction dissolves rapidly. 

We aim to analyze the relation between field stars and 
star clusters in different environments and address the ques- 
tion of whether or not there is a constant cluster formation 
"efficiency". To this end, we use Hubble Space Telescope 
(HST) images of a set of five galaxies (NGC 4395, NGC 
1313, NGC 45, NGC 5236, and NGC 7793), which are 
nearby, face-on, spirals that differ in their morphologies, 
star, and cluster formation histo ries. The cluster sy stems 
of these galaxies were studied by iMora et all (|2009D using 
the same Hubble images analyzed in this work. Mora et al. 
observed significant variations in the cluster age distribu- 
tions of these galaxies. The galaxies are sufficiently nearby 
(^ 4 Mpc) for the brighter field stars to be well resolved in 
HST images, so that (recent) field star formation histories 
can be constrained by means of the synthetic CMD method. 
We can therefore take advantage of the superb spatial res- 
olution of HST images to study field stars and clusters si- 
multaneously within specific regions of these galaxies. 

As a pilot work, this paper is devoted to presenting and 
testing all the analysis procedures, such as detection of stars 
and cluster candidates, completeness tests and photometry, 
and derivation of field star and cluster age distributions. We 
discuss our implementation of the synthetic CMD method 
as an IDL program and carry out tests of this program. 
To test our procedures we use the galaxy NGC 4395. The 
methods described in this paper will be used in our study 
of the rest of the galaxy sample (Silva- Villa et al. 2010, in 
prep.). 

This article has the following structure. In Sect. 2, we 
provide general information about NGC 4395. We present 
the observations, data reduction, and photometry, where 



we differentiate the point sources (star candidates) from 
the extended objects (cluster candidates) in Sects. 3 and 4. 
Section 5 is devoted to the analysis of the star and cluster 
properties. Finally, in Sect. 6 we present the discussion and 
conclusions. 



2. NGC 4395 

According to the NASA/IPAC Extragalactic Database 
(NED), NGC 4395 is a late- type spiral classified as type 
SA(s)m. It harbours the closest and least luminous known 
example of a Seyfert 1 nucleus (iFilippenko fc Sargent! 
119891 ). Following iLarsen fc Richtle^ (|1999D . we adopt a dis- 
tance modulus of (m - M) = 28.1 {D « 4.2 Mpc) and 
an absolute magnitude Mb — —17.47, intermediate be- 
tween the Small and Large Magellanic Cloud. We assume a 
Galactic foreground e xtinction for NGC 4395 of Ab = 0.074 
(ISchlegel et al.lll998h . 




Fig. 1. Second Palomar Sky Survey images of NGC 4395 
combined using Aladin software. HST/ ACS (white lines) 
and HST/WFPC2 (yellow fines) fields covered by our ob- 
servations are indicated. 

The cluster system of N GC 4395 was first studied 
by iLarsen fc RichtleJ (|1999|) . using ground-based imag- 
ing. Using multiband (UBVRIHa) photometry, Larsen fc 
Richtler identified 2 young clusters in this galaxy, although 
their observations were limited to objects with My < —8.5. 
In their work, NGC 4395 is part of a sample of 21 galaxies. 
Compared to the remaining sample, NGC 4395 exhibits 
an exceptionally small number of star clusters and a low 
specific luminosity (rati o of cluster- to total g alaxy light; 
Larsen fc Richtler) [2000l ). Using HST images, iMora et all 
(2009) could detect clusters based on their sizes and found 
a total of 44 clusters in NGC 4395 to a magnitude limit of 
Mb = —3.3. Compared to the other 4 galaxies in their sam- 
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Table 1. Journal of HST/ACS and HST/WFPC2 observations for both fields (Fl and F2) in NGC 4395. 



Proposal ID Date Filter 


Total exp. time (s) 


NGC 4395 (J2000.0) a : 12''26™00'' 5 : +33°3l'04" 
9774 2004 Jun 12 F435W (B) 
9774 2004 Jun 12 F555W (V) 
9774 2004 Junl2 F814W (I) 


680 
680 
430 


NGC 4395 (J2000.0) a : 12''26'"00".50 5 : +33°30' 58" .3 
9774 2004 Jun 12 F336W (U) 


2400 


NGC 4395 (J2000.0) a : 12''25'"45'.20 5 : +33°34'28" 
9774 2004 Jun 18 F435W (B) 
9774 2004 Jun 18 F555W (V) 
9774 2004 Jun 18 F814W (I) 


680 
680 
430 


NGC 4395 (J2000.0) a : 12''25'"42".72 5 : +33°34'22".6 
9774 2004 Jun 18 F336W (U) 


2400 



pie, Mora et al. reached the same conclusions as Larsen & 
Richtler, i.e., NGC 4395 has a small number of star clusters. 
Mora et al. estimated the ages and masses of the clusters 
detected, showing that the star cluster system contains no 
objects with ages older than 10^ yr and includes clusters 
with masses ranging from 200 to ~ 3 x 10^ M©. 

Until now, no study of resolve d field stars has be e n per - 
formed in this galaxy. However, 'Lars en fc Richtleil (|2000[ ) 
found that NGC 4395 has one of the lowest area-normalised 
star formation rates of all objects in their sample. 



3. Observation and data reduction 

Images of NGC 4395 were taken using the Wide Field 
Channel on the Advanced Camera for Surveys (ACS/WFC) 
and the Wide Field Planetary Camera 2 (WFPC2), both 
on board the Hubble Space Telescope (HST). The resolu- 
tion of the detectors are 0'.'05, 0'.'046, and O'.'l per pixel for 
ACS/WFC, WFPC2/PC, and WFPC2/WF, respectively. 
At the distance of NGC 4395, 0'.'05 corresponds to a linear 
scale of ^ 1 pc. 

Two different fields of the galaxy were observed, cover- 
ing the two spiral arms (see Fig.[T]). The images were taken 
using the fihers F336W (- U) using WFPC2 and F435W 
{~ B), F555W (- V) and F814W (~ I) using ACS, for each 
field. Each exposure was divided into two sub-exposures to 
eliminate cosmic-ray hits. For the ACS images, these sub- 
exposures were also dithered to allow removal of cosmetic 
defects. Table [1] summarizes the observations. 

The data were processed with the standard STScI 
pipeline. The raw ACS images were d rizzled using the mul- 
tidrizzle tas kjlKoekemoer et aI]|2002D in the STSDAS pack- 
age of IRAfQ. The default parameters were used, but au- 
tomatic sky subtraction was disabled. The WFPC2 images 
were combined and corrected for cosmic rays using the crrej 
task with the default parameters. 

3.1. Object detection 

At the distance of NGC 4395, the spatial resolution of 
the HST images allows us to distinguish field stars (point 



^ IRAF is distributed by the National Optical Astronomical 
Observatory (NOAO), which is operated by the Association of 
Universities for Research in Astronomy, Inc, under cooperative 
agreement with the National Science Foundation 



sources with typical FWHM ^ 2 pixels) and star clus- 
ters (extended sources with typically FWHM ~ 3 pixels 
or greater) in the images. We analyzed the two separately 
and applied different detection criteria optimised for each 
type of object: 

— Field star candidates (point sources) 

To detect field stars, we created an averaged image (us- 
ing the bands B, V, and I) for each field. We ran the 
daofind task in IRAF for the detection of stars using a 
4ct detection threshold and a background standard de- 
viation of ~ 0.02 in units of counts per second. 

— Star cluster candidates (extended sources) 

Cluster detection was also performed on the aver- 
aged image. The ob ject detection was perfor med using 
SExtractor V2.5.0 |Bertin &: ArnoutsI Il996h . The pa- 
rameters used as input for this program were 6 con- 
nected pixels, all of them with lOcr over the back- 
ground, to remove point- and spurious sources as much 
as possible. From the output file of SExtractor, we 
kept the coordinates of the objects detected and the 
FWHM calculated by this program. The coordinates 
found wit h SExtractor were then passed to ishape in 
BAOLab (jLarsei]|ll999r) . Ishape models star clusters as 
analvtical iKind ( |l962f) profiles (with concentration pa- 
rameter tfidai/fcore = 30) and takes the instrument's 
PSF, created over the average image during the pho- 
tometry, into account (see Sect. 4 for details on the cre- 
ation of the PSF). By minimization of a x^-like func- 
tion, ishape calculates the best-fit cluster coordinates, 
size, i.e., FWHM, signal-to- noise ratio, and the of 
the best fit. From the output of ishape we saved the 
coordinate of the objects, the FWHM, the of the 
fit, and the signal-to-noise ratio (calculated within the 
fitting radius of 4 pixels) . 

For each WFPC2 chip, between five and seven common 
stars were visually selected and used to convert the coor- 
dinate list from ACS to WFPC2 frames using the task 
geomap in IRAF. The transformations had an rms of 
~ 0.15 pixels. 



4. Pliotometry 

The photometry of NGC 4395 was performed following 
standard aperture and PSF fitting photometry procedures 
as described in the following. 
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Fig. 2. Magnitude errors for the stars detected using PSF 
fitting photometry over the two fields for the bands B, V, 
and I 

4. 1 . Field stars 

Because of the crowding in our fields, we performed PSF 
fitting photometry to study the field stars. A set of bona fide 
stars were selected by eye to construct the PSF (for each 
band a different group of stars were used because the same 
stars might have different brightnesses in different bands). 
The PSF photometry was performed with DAOPHOT in 
IRAF. We selected the PSF stars by measuring the FWHM 
(using imexamine) and selecting point sources smaller than 
FWHM fa 2.2 pixels. As far as possible, we tried to include 
isolated stars distributed over the whole image. 

The raw magnitudes were converted to the Vega mag- 
nitude system using the HST zero-points taken from HST 
webpagef0 after applying aperture corrections to a nominal 
O'.'S aperture (see sect. 4.4 for a more detailed description 
of how aperture corrections were determined). The zero- 
points used were ZPb = 25.767, ZPy =25.727, and ZPi = 
25.520 magnitudes. 

Figure m shows the errors in our PSF photometry {la 
error) versus magnitude. The errors increase strongly below 
magnitudes of ~ 26 in each band, corresponding to absolute 
limits of ^ —2 mag at the distance of NGC 4395. 

A total of ~30 000 stars were found in each field. A 
combined Hess diagram for both fields is shown in Fig. [3] 
(a Hess diagram plots the relative frequency of stars at dif- 
ferent color- magnitude positions). Various phases of stellar 
evolution can be recognised in the Hess diagram: 

1. Main sequence and possible blue He-core burning stars 
at 1/ - / ~ and -2 < F < -8; 

2. Red He core burning stars at 1.2 < V — I < 2.5 and 
-2.5 <V < -6.5; 

3. RGB/AGB stars at 1 < F - / < 3 and -0.5 <V< 
-2.5. 

Overplotted in Fig.[51are theoret ical isochroncs from the 
Padova group (jMarieo et al.l 120081 ) for five different ages 
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Fig. 3. Hess diagram for the field stars in both fields. The 
dashed white line represents the 50% completeness for the 
first field. Red lines are Padova 2008 theoretical isochrones 
for the ages 10'^,10^■^10^10*■^ and 10^ yr using LMC 
metallicity. 



using metallicity Z = 0.008 (red lines). The O abundance 
of Hli regions in NGC 4395 was measured by iRov et all 
(1996) to_bel2+]ogO/H = 8.33 ± 0.25, or about 1/3 solar 
(Grcvcss e fc Sauvall 19981 ). suggesting an overall metaUicity 
similar to that of the LMC. This is consistent with our 
analysis in Sect. 5, which shows that isochrones of LMC 
metallicity reproduce our data more closely. The white line 
indicates the 50% completeness limit (see Sect. 4.3). 

4.2. Star clusters 

Our detection criteria is met by 16 463 objects, i.e., 6 con- 
nected pixels with IOct above the local background level. 
For the clusters, we carried out aperture photometry. An 
aperture radius of 6 pixels was used for the ACS images. At 
the distance of NGC 4395, 1 ACS/WFC pixel corresponds 
to ~ 1 pc, hence the chosen source aperture contains about 
2 half-light radii for a typical star cluster. Our sky annulus 
had an inner radius of 8 pixels and a width of 5 pixels. For 
the WFPC2 images, apertures covering the same area were 
used (source aperture = 3 pixel radius, sky annulus = 4 
pixel inner radius and 2.5 pixels width). 

Of the 16 463 objects detected with SExtractor/ is/iape, 
a total of 4 472 candidates have measured four band pho- 
tometry. 

ACS magnitudes were converted to vega magnitude sys- 
tem using the same tables as for the field stars. WFPC2 
magnitudes were converted to the Vega magnitude system 
using the zero-points taken from the webpage^ of HST. 



^ http:/ /www.stsci.edu/hst/acs/analysis/zeropoints/^^tablestart ^ www.stsci.edu/instruments/wfpc2/Wfpc2_ch527^1933986 
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Fig. 4. Histograms of the FWHM of the objects detected 
with SExtractor (top panel) and ishape (bottom panel). 
The vertical dotted lines represent the limits used to se- 
lect the extended objects [FWHMsExtr actor > 2.7 and 
FWHNUshape > 0.7 pixels) 

Charge transfer efficiency ( CTE) c orrections were applied 
following the equations from iDolphi n (200(Jj3. 

To select star cluster candidates, we used 3 criteria: 

— Size: 

We measured the sizes of the objects. Figure |4] shows 
the FWHM distributions for SExtractor and ishape. 
The SExtractor histogram peaks at ~ 2.2 pixels, cor- 
responding to the stellar PSF. The ishape histogram 
peaks at pixels, as ishape takes the PSF directly into 
account. Based on the FWHM distributions, we decided 
to use a criteria of FW HMsExtr actor > 2.7 pixels and 
hape ^ 0-7 pixels as a first selection of ex- 
tended sources. These two limits correspond to a phys- 
ical cluster half-light radius of 1 pc or greater at the 
distance of NGC 4395. Many of the candidate clusters 
we detect are low-mass and have irregular profiles often 
dominated by a few stars, so we chose to rely on both 
SExtractor and Ishape size measurements to achieve a 
more robust rejection of unresolved sources. 

— Magnitude: 

ora et all (|2009l ) ran a completeness analysis of arti- 
ficial star clusters of different FWHMs in five nearby 
galaxies, including NGC 4395. The results presented in 
their work establish a 50% limit between 25 < tob < 
26 magnitudes f or objects w i th FWHM = [0.1 (point 
sources), 1.8] fsee lMora et"all (|2007L l2009t) and Sect. 4.3 
for more details). 

By performing four band photometry, we set a mag- 
nitude cut-off at mv < 23, which represents objects 
brighter than My ^ — 5 at the distance of the galaxy. 

^ L ast update May 12,2008; 

|http : //purcell . as . arizona. edu/wfpc2_calib/| 




Fig. 5. Example stamps of the clusters that satisfy the cri- 
teria stated in this work. Each cluster image has a dimen- 
sion of 100 X 100 pixels (- 100 x 100 pc). 

This limit is brighter than that found by iMora et al.l 
( 2007) by 2 magnitudes. 

The main differen ce between the selection criteria of 
iMora et al.l (|2007D and this work is in the size lim- 
its of ishape, of half of a pixel (Mora et al. used 
FWHMishapc > 0.2), hence we expect our data not to 
be significantly affected by incompleteness to our mag- 
nitude cutoff. However, even at our limit of My = — 5, 
there is a risk that a few stars may dominate the light 
originating in a cluster, making it difficult to differenti- 
ate a real cluster from a couple of stars that, by chance, 
could be in the same line of sight. 
- Color: 

Without taking into account any significant reddening, 
all clusters, including globulars, w i ll hav e colors bluer 
than V-I ^ 1.5, e.g., Forb es et all (|l997^ ■ lLarsen et aP 
(2001). We therefore make a color cut dXV — I = 1.5. 

We wish to emphasise that there is probably no unique 
combination of criteria that will lead to the detection of all 
bona-fide clusters in the image and at the same time pro- 
duce no false detections. At low masses and young ages in 
particular, the light profiles may be dominated by individ- 
ual bright stars. Out of 4 472 candidate objects, a total of 
22 objects fulfill the three criteria stated above and will be 
considerate in the remain of this study to be star clusters. 
These 22 clusters were visually inspected in our images to 
determine whether they resemble a cluster. An example of 
the clusters detected is presented in Fig. [S] We investigated 
the large number of rejected objects and found that they 
were rejected for many reasons: (1) the area covered by the 
two detectors differs by a factor of ~ 2; (2) many of the ob- 
jects are too faint in the U band; (3) the magnitude cut-off 
removes many objects, e.g., from a total of 4 475 detected 
objects, only 1 351 satisfy the magnitude limit, leading to 
the loss of ~ 70% detections; and (4) a large fraction of the 
objects have FWHMs smaller than our limits. 

For the extended objects, a two-color diagram, based 
on the photometry in all of the available 4 passbands, is 
presented in Fig. |6l Clusters in both fields are depicted 
with their respective photometric errors and corrected for 
foreground extinction [Ab = 0.074). Using GALEV models 
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Fig. 6. Two-color diagram for the star clusters detected in 
both fields, x represents the first field and A the second 
field. The lines represent the theoretical GALEV track for a 
cluster with LMC (dashed line) or solar (dash-dotted line) 
metallicity, adopting Padova isochrones, a Salpeter IMF, 
and ages between 6.6 < log{T/yr) < 10.1 



( Anders fc Fritze-v . AlvenslebenI I2003D , Padova isochrones 
(jBertelh et al.lll994[ ). a Salpeter's IMF (iSalpeteijll955l ). and 
LMC (dashed line) or solar (dash-dotted line) metallici- 
ties, we overplotted the track that a cluster follows from 
ages between 6.6 < Log{T) < 10.2 yr (each A log(r) — 0.5 
dex age in log units being indicated) . We observe consider- 
able scatter in the two-color diagram. For these relatively 
low-mass clusters, the discreteness of the initial mass func- 
tion will be important and might contribute significantly to 
the scatter ( Girardi et al.lll995l : [Cervino fc Luridianall2006t 
iMai'z Apellan iz 2009). A more detailed study of stochastic 
effects will be presented in a forthcoming paper (Silva- Villa 
et al. 2010, in prep.) This scatter was previously observed 
bv iMora et al.l (|2009) (see their Fig. 6). Considering this 
scatter, it is clear that ages derived by a comparing ob- 
served and model colors should be treated with some cau- 
tion. 



4.3. Completeness test 

To determine the completeness limits of our field-star pho- 
tometry, we generated synthetic images with artificial stars. 
From 20 to 28 magnitudes, in steps of 0.5 mag, five images 
per step were created and analyzed with exactly the same 
parameters used in the original photometry. For each com- 
bination of color and magnitude, 528 artificial stars were 
added using 4 different regions over the image to cover 
crowded and uncrowded areas. These test images were cre- 
ated using the task mksynth in BAOLab (Larsen 199^ and 
using the original PSF images created during the photom- 
etry procedures (see Sect. 4.1). The separation between 




Fig. 7. Subsection of the original image (right) and com- 
pleteness image (left). 



two consecutive stars was 100 pixels (without any sub-pixel 
variations), avoiding possible overlaping among the stars. 
The resulting test images were added to the science images 
using the task imarith in IRAF. As an example, a subsec- 
tion of an image, both test and science, is presented in Fig. 
[71 where the fake stars have magnitudes mv = 21. 

To performe a realistic completeness analysis it is 
in principle necessary to sample the three-dimensional 
{B,V,I) color space. However, since different colors are 
tightly correlated with each other, the problem can be re- 
duced to a two-dimensional one. Figure [S] was generated 
using Padova 2008 isochrones (Marigo ct al. 20Q8|), assum- 
ing solar metallicity, in the color range from -1 to 2, for 
B — V and V — I and shows that there is a nearly 1:1 
relation between these two colors (the red line in Fig. |5] 
represents a 1:1 relation, but not an accurate tof the data). 
We created the images for the completeness test described 
above using this approximation between the three bands, 
i.e., if a i? image has stars with tob = 21 and B — V = 1, 
then the V and / images will have stars of my = 20 and 
mi = 19, respectively, allowing us to perform a study of 
magnitudes and color variations over the images. 

Based on these tests, we found an average of 50% com- 
pleteness limits for the whole color range at = 26.69, 
mv — 26.55, and m/ = 26.42 for the first field and 
tob = 26.71, my = 26.49, and mj = 26.39 for the second 
field. Figure [3] shows the completeness diagram obtained 
from this analysis for the first field (for the second field, 
the figure is similar) and each magnitude. 

The color-dependent 50% completeness limit is shown 
as a white dashed line in the Hess diagram (see Fig. [3]). 
Since the completeness functions are very similar for the 
two fields, we can combine the photometry for both fields 
and use just one set of completeness tests in the following 
analysis. 

We did not performe a completenes s test for star clus- 
ters b ut refer to the tests performed bv lMora et al.l ()2007l 
I2OO9D . who used the same data and very similar cluster 
detection procedures. These authors created artificial star 
clusters using different FWHMs from O.lpixels (stars) to 
1.8 pixels and a range of magnitudes from 16 to 26 for 
three square grids in different positions over the image, 
trying to cover crowded and non-crow ded ar e as. T o cre- 
ate the fake clusters, they assumed a iKind (|1962D pro- 
file with Ttidai/Tcore = 30. Using the mkcmppsf task on 
BAOLab, fake extended objects were added to an empty 
image and then added to the science im age to pe r forme 
meas urements. The selection criteria in Mora et al.l (|2007L 
I2009D is a FWHM > (2.7, 0.2) pixels for SExtractor and 
ishape, respectively. For high background levels. Mora at al. 
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Table 2. For point source. B, V, and I aperture corrections 
to a nominal aperture of 0."5, estimated in this study. 



Fig. 8. Two color diagram for theoretical values of B, V 
and I bands from Padova 2008 isochrones adopting LMC 
metallicity for ages between 10®'^ to 10^" yr. The red line 
represents a 1 : 1 relation between the colors B-V and V-I. 



found a shallower detection limit; nevertheless, all the lim- 
its correspond a 50% completeness limit between « 25 
and TUB ~ 26, 2-3 magnitudes fainter than the cut at V=23 
that we apply for the selection of cluster candidates. 
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Fig. 9. Completeness diagrams for the bands B,V, and I 
studied over the first field. Vertical lines represent the 50% 
completeness. 
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4.4. Aperture corrections 

The aperture corrections were determined separately for 
field stars and star clusters. The aperture corrections for 
the field stars were derived following standard procedures, 
w hile for extended ob jects we adopted the relations found 
bv lMora et all (|2009D . 



Field stars: 

By "aperture corrections" , we here mean corrections 
from the PSF-fitted instrumental magnitudes to aper- 
ture photometry for nominal radii of O'.'S. These were 
measured using a set of isolated visually selected stars 
across the images. F rom 0'.'5 to infinity, we applied the 
ISirianni et al.l (|2005l) values. The corrections obtained 
are of the order ~ 0.1 mag (see table [2|). 
Star cluster s: 

[M ora et all ()2009D estimated a relation between the 
aperture corrections and the sizes (FWHM) of star 
clusters using the same data s et used in this pa per. 
Photometric parameters in both iMora et al.l (|2009l and 
ours) work are the same, allowing us to assume the re- 
lations found in their work and apply these aperture 
corr ections (size-depend ent) to our data, following eq. 
1 in ora et all (j2009[ ). This set of equations is also 
band-dependent, although we used the sizes of the ob- 
jects measured on an average image . 
The aperture corrections by iMora et al.l (l2009f l corre- 



spond to a nominal aperture of 1'.'45. From this nomi- 
nal apert ure to in finity, we adopted the values presented 
by Sirian nTet al.l (f2005), although these corrections are 
^ 0.03 magnitudes for the bands B, V, and I (within 
1'.'5 about 97% of the total energy is encircled). 



5. Star formation histories, ages, and masses 

Our main goal in this paper is to compare the field stars 
with the star cluster populations in NGC 4395. In the fol- 
lowing section, we describe how we derive the star forma- 
tion histories (SFHs) of the field stars and the ages and 
masses of the clusters. We then proceed to compare the 
two. 



5.1. Deriving star formation histories: approacli and its 
testing 

To estimate the SFH of the field stars, we have used the 
synthetic CMD method. The synthetic CMD method con- 
sists of creating of an artificial photometric distribution of 
stars taking into account photometric errors, distance mod- 
ulus, IMF, binarity, and interstellar extinction. The relative 
weights of the stellar isochrones used in generating the syn- 
thetic CMD are adjusted until the closest possible match of 
the synthetic CMD to the observed one is obtained. These 
weights can then be translated to a SFH. 
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Fig. 10. SFH test using the bands B, V, and I. Artificial populations were created using Padova and Geneva isochrones 
and then passed to the program. The procedure was performed 10 times and the mean SFH (dashed line) and the 
respective standard deviation is shown, relative to the input assumed SFH of OT Msti„yr~^ (straight line). Each row 
represents a different binning used for the study, which decreases in size from top to bottom. 



F ollowing the wo r k oflTbsi et al.l (|l99lf ). lDolDhh] (|1997D . 
and lAparicio et al.l ()1997[ ). many authors have imple- 
mented this technique. Some examples are STARFISH by 
iHarris fc Zarit skv (2002), I A C- Star by lAparicio fc GallartI 
(I2004D. lAC^op by Aparicio fc Hidalgol (|2009D . among oth- 



ers. iSkillman fc Gallarti (|2002D presented the results of a 
workshop (The Coimbra Experiment) devoted to compar- 
ing different methods and interpretations of the results from 
various groups using a homogeneous data set and physical 
inputs. The results showed a large scatter at young ages, 
but good agreement at old ages. 

We implemented the synthetic CMD technique our- 
selves as an IDL program. We chose to develop our own 
code to include any features that we might wish (or need), 
and there is no clear preference in the literature for any 
of the other codes. The IDL code used in this work was 
based on o ur previous experie nce with the synthetic CMD 
technique (|Larsen et al.l |2003) . The program searches for 
the best-fit synthetic Hess diagram using a maximum- 



likelihood technique. We assume that the likelihood of ob- 
serving ki stars in the ith pixel of the Hess diagram is given 
by the Poisson distribution 'P{ki] iii), where fii is the num- 
ber of stars expected in this pixel according to a given 
model Hess diagram. The program then searches for the 
linear combination of input isochrones that maximises the 
total likelihood C over the entire Hess diagram 



log/: ^^hgVih-^n,). 



(1) 



In practice, we assign a low (non-zero) probability of having 
a star even in pixels where fii = 0. This is essentially equiva- 
lent to allowing a constant "background level" of stars that 
are not fitted by any isochrone. This is necessary since any 
star in a pixel with fii = would otherwise immediately 
drive the total likelihood to zero. Since no combination of 
existing isochrones provides a perfect match to the observed 
CMDs, and one may in any case always expect some con- 
tamination of the CMD by foreground stars, background 



E. Silva- Villa and S. S. Larsen: The star cluster - field star connection in nearby spiral galaxies 



9 




-0.95 



° -1.05 



-0.95 



-1 ,05 



7.0 



HH- 



B-V 



III 



SFR=0 
CMPL=-0, 
CMPL=0 
CMPL=0, 
CMPL=1 



.5 8.0 

bg(age) [Yr] 



9.0 



Fie. 11. Hess diaeram for the 4 possible completeness limit tt i r, rn , r , i cmrT • j-n j. i j. i- 

, . ppS ^ ^ Fig. 12. lest or the bhti usmg dirterent completeness lim- 



studied in Fig. [H 



its. 



galaxies, etc., this would be overly restrictive. We also ex- 
perimented with other criteria, e.g. minimum r.m.s. differ- 
ence between observed and model Hess diagrams, but tests 
similar to those described in Sect. 5.2 below indicated that 
the maximum-likelihood method gave the most reliable re- 
sults. 

To create the synthetic CMD, the program uses theo - 
retical isoc hrones either from Padova (|Marigo et al.|[2008h 
or Geneva (jLeieune fc SchaereH[200l or any other set of 
isochrones, as long as they tabulate the relevant color ver- 
sus mass. If the weight of each isochrone were fitted in- 
dependently, this would lead to a very large number of 
free parameters, due to the small difference in age between 
individual isochrones, e.g. Alogr = 0.05 for the Padova 
isochrones. We therefore group the isochrones together in 
age bins that can be defined by the user, typically using a 
bin size of Alogr = 0.15 dex. 

The program allows us to assign different weights to 
different rectangular regions of the CMD, selected by the 
user. This is useful because some phases of stellar evolution 
are more uncertain than others, e.g. the blue loop stars, and 
thus should carry less weight in the overall determination of 
the SFH. The use of completeness limits (if known) can be 
applied before the creation of the synthetic models, inter- 
polating the completeness limits over the whole area of the 
CMD. The option of applying the completeness correction 
after the creation of the synthetic model is simply a multi- 
plication of each pixel in the synthetic CMD with a number 
between and 1. We also include a simplified treatment of 
binaries, in which binary evolution is ignored but the ef- 
fect of unresolved binaries on the CMD are modeled. The 
program currently allows three different assumptions about 
the mass ratios in binary systems, namely a delta function, 
an IMF sampling and a flat distribution. The metallicity, (a 
range of) extinction values, and distance modulus must also 
be specified. We allow the extinction to be age- dependent. 



by providing a list of ages and a range of extinction values 
for each age. The program then interpolates in this list for 
each isochrone. 

Isochrones are populated by assuming a lSalpeteil ()1955() 
IMF with SFHs normalised to a mass range specify by the 
user (M = [0.15,100] Mq is the default range). Since we 
use Hess diagrams, each point in the CMD is a density func- 
tion. To construct them, different kernels, namely square, 
Gaussian, disc, or delta function, can be used. The kernels 
are of adjustable resolution and dimension. In the test pre- 
sented below, we adopted a resolution for the Hess diagrams 
of 100 X 100 pixels and used a delta function kernel. Each 
isochrone is broadened by the assumed photometric errors 
and binarity, and shifted and broadened by the (range of) 
specified extinction values. To ensure a smooth coverage 
of the CMD, the program interpolates the isochrones by a 
factor of 10. After creating of the Hess diagram for each 
individual isochrone, the program linearly combines them 
and a synthetic CMD Hess diagram is obtained. 

We performed a number of tests to check the internal 
consistency of the program, as well as the sensitivity of the 
derived SFHs to the various parameters involved. 

5.2. Internal consistency check of the code 

We created artificial stellar populations, passed them to 
the program, and then checked how the output compared 
to the input. The populations were constructed assuming a 
constant SFR of 0.1 Moyr"^, Padova or Geneva theoreti- 
cal isochrones, a Salpeter's IMF, and solar metallicity. The 
mass of each star was randomly chosen using Salpeter's pre- 
scription for masses between 1 and 100 Mq (low-mass stars 
(< 1 Mq) are too faint to appear in the Hess diagram, 
especially after applying a completeness limit). Because 
we restricted the mass range in our analysis, we then ex- 
trapolated to the default mass range used by the code 
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Fig. 13. Padova theoretical isochrones were used to test the SFH program developed using the photometric data obtained 
in Sect. 4. Panel a): Test for different metallicities (solar, LMC, SMC) without binarity or extinction assumed. Panel b): 
Test for different extinctions assuming LMC metallicity and no binarity. Panel c): Test for different binarity assuming 
LMC metallicity and no extinction. Panel d): Test for different boxes used over the fitting data assuming LMC metallicity, 
no extinction, and no binarity. 



{M — [0.15,100] M0). The age of each star was assigned 
randomly (from a uniform distribution) for ages between 4 
Myr to 1 Gyr. Having the mass and the age of each star, 
magnitudes were obtained by interpolating the theoretical 
isochrones. The magnitude limit used in these tests was 
My = 0. We did not include binaries or extinction in this 
initial internal consistency check of the code. 

Figure [TUl presents the results of this analysis. In all the 
panels, the input SFH is shown as a horizontal straight line 
and the average reconstructed SFH of the output after 10 
different runs as dashed lines. Each run had a new popula- 
tion, i.e. different random realization, for the bands B, V, 
and L The error bars illustrate the standard deviation in 
the reconstructed SFHs. We studied the variations in these 
results based on the different assumptions of binning. Three 
different bins were used in this test corresponding to each 
row of Fig. [TU] 



We see that the reconstructed star formation histories 
generally agree fairly well with the input values regardless 
of the binning, color, or isochrones used. It is observed that 
the error bars are larger for smaller bins, although good 
agreement with the input SFH is maintained. The maxi- 
mum difference is Alog{S FR) ^ 0.05 dex, and generally 
far less. 

We performed a test to see how completeness affects 
our SFH estimation. Figure [TT] shows the Hess diagram of 
a stellar population with ages between 10^ to 10^ yr using 
the mass range from above. Overplotted are 3 Padova 2008 
isochrones for ages 10^, 10^, and 10® yr and four complete- 
ness limits at My = —0.5,0.0,0.5, and 1.0 magnitudes. 
Figure [T^ shows that the SFHs can be recovered reliably to 
progressively older ages as the completeness limit becomes 
fainter. In particular, the "burst" at 10^'^ years disappears 
when fainter stars are included in the fit. 
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Fig. 14. Comparison of the fit for different nietallicities (so- 
lar, LMC, and SMC). 
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Fig. 15. Bootstrapping study. The errors retrieved from 
this test are the random errors of our program. No strong 
variations caused by random errors affect our results. 



ICignoni fc Tosil (|2010f ) created artificial CMDs to simu- 
late a stellar population with a constant SFR between — 13 
Gyr and used four different completeness limits. They con- 
cluded that to safely reconstruct a SFH over a full Hubble 
time from a CMD, we need to resolve all stars down to 
My = 4.5. However, this completeness limit can only be 
reached for galaxies in the Local Group. 

We conclude that our program produces consistent re- 
sults at the level of 0.05 dex, the age range over which this 
holds being restricted mainly by the completeness limits. 
This internal consistency of the code does not, however, 
imply that real SFHs can be recovered with similar accu- 
racy. 

5.3. Tests of sensitivity to assumptions about parameters 

We now concentrate on testing the effects of four param- 
eters: metallicity, extinction, binarity, and the fitting area 
used (rectangular boxes). To perform this test, we executed 
the code using the scientific data described in Sect. 4, using 
different assumptions about these parameters and a mass 
range M = [0.10, 100] Mq (which is the same mass range 
used in Sect. 5.4). The results are shown in Fig. [T21 Each 
panel represents one of the tested variables. 

In panel a), we show the SFHs derived for three assump- 
tions about metallicity: solar (Z=0.02), LMC (Z==0.008), 
and SMC (Z=0.004). In this panel, we ignore binaries, in- 
clude no additional extinction and used the fitting area 
called "boxes2" described below. We observed that the SFH 
does not change very much from LMC to solar, but at SMC 
metallicity there is a much stronger increase from ^ 100 
Myr ago to the present. By comparing the observed and 
model Hess diagrams, it is clear that the solar metallicity 
isochrones generally provide a poor fit, especially for the red 
and blue supergiants that appear much too cool compared 



to the observations, while LMC and SMC metallicities are 
very similar (see Fig. [14]). 

Panel h ) presents the behavior of the SFHs for three dif- 
ferent assumptions about the total extinction (foreground 
plus internal): Af, = [0.1,0.15,0.2]. We fixed the metallicity 
at Z = 0.008 (LMC), again ignoring binaries, used Padova 
isochrones and used the fitting area "boxes2". In general, 
the SFHs obtained using different extinctions do not vary 
strongly for these relatively modest variations in the ex- 
tinction. Greater extinction variations would produce larger 
shifts towards the red in the model Hess diagrams and are 
inconsistent with the data. 

Panel c) shows the dependence on the assumptions for 
binary fraction (/) and mass ratio (g) of values (1.) / = 0, 
q = 0; (2.) / = 0.5, q = [0.1,0.9]; (3.) / = 1, q = 1. The 
panel also indicates how this dependence affects the esti- 
mation of the SFHs. We used Padova isochrones, fixed the 
metallicity at Z = 0.008, ignored the effects of additional 
extinction, and used the fitting area named "boxes2" de- 
scribed below. Binarity can have an effect on the derived 
SFHs with minor effects for r < 100 Myr and a shift to- 
wards somewhat lower overall SFR for the extreme case 
/ and q equal 1. However, the binarity assumptions made 
here are not completely realistic. The first and third as- 
sumptions are extreme cases with no binaries at all and 
a galaxy where all the stars have a companion of exactly 
the same mass, respectively. The second assumption is an 
intermediate case where there is a continuous (flat) range 
of q values and a star has ^ 50% probability of being in a 
binary system. However, in a more realistic case binary evo- 
lution should be taken into account as well. This is however 
beyond the scope of the present work. 

In the last panel, panel dj, we show the SFH esti- 
mates obtained using two different sets of boxes in the 
fltting. As overplotted in Fig. [T71 where we plot the ob- 
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served Hess diagram and the regions fitted, the left box 
(-0.5 < V - I < -0.05 and -8.0 < V < -2.5) used 
contains the main sequence stars, the blue He core burn- 
ing phases being outside of the fit, while the right box 
(1.0 <V -I <2.5 and -8.0 <V < -2.0) contains the red 
He core burning stars and possible RGBs and AGBs. This 
set of boxes is called " boxes2" . The second set of boxes cov- 
ers the same area for the evolved stars (right box), while 
the blue He core burning stars are also included in the fit- 
ting box (-0.5 <V - I < 0.7), which we label "boxesl" 
in panel d). In general, we found that, the SFH result from 
"boxes2" shows a lower value that for "boxesl" . We studied 
these two possible sets of boxes because the fitted Hess dia- 
grams retrieved using these two boxes indicated that the fit 
for the blue He core burning stars is not in good agreement 
with observations. The variations after ~ t ^ lO*'^ yr are 
unlikely to be real, but are probably due to incompleteness 
as discussed above. 

We created a fake population using Padova isochrones 
and assuming solar metallicity, no extinction, and no bina- 
rity; we then analyze this with our program, using Geneva 
isochrones to see whether the assumed SFH could be recov- 
ered. Figure [16] shows the result of this test. There are sig- 
nificant differences between the input and recovered SFH, 
deviations being as large as ~ 0.4 dex at young ages (< 10 
Myr) and differences being at the level of 0.1-0.2 dex at 
older ages. The error bars presented in Fig. \W\ were deter- 
mined using the same method used in the previous section, 
i.e., are equivalent to the standard deviation of the recon- 
structed SFHs after 10 different runs with different random 
realizations. 

We finally test how random errors (due to the finite 
number of stars in the Hess diagram) affect the SFHs, us- 
ing a bootstrapping method whereby the input data, i.e., 
the photometric data found in this work, are randomly re- 
sampled many times. To perform this test, we used the ob- 
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Fig. 16. Test using Padova isochrones to create the popu- 
lation and Geneva isochrones to analyze them. 
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Fig. 17. Gomparison between the observed and fitted Hess 
diagrams for the field stars studied in this work. As age 
indicators, three isochrones are overplotted at the ages of 
10^, 10^, and 10^ yr. Boxes used to fit the data are repre- 
sented for both Hess diagrams, observed and fitted. 



servations obtained in this work (see Sect. 4), Padova 2008 
isochrones, LMC metallicity, no additional extinction, and 
no binarity. Disabling binarity will not influence our results 
dramatically and will save us computational time. We ran 
the program 100 times, redistributing randomly the pho- 
tometric data in each run. From Fig. I15[ we conclude that 
Poisson variations do not significantly affect our results. 

To summarize, we demonstrated that the main uncer- 
tainties in our derived SFHs are systematic, depending on 
assumptions about binarity, the choice of isochrones, and 
metallicity. However, with the exception of the use of so- 
lar metallicity isochrones (which provide a poor fit to the 
data) , the overall mean SFRs are fairly consistent with our 
data up to ages of a few hundred Myr. 

5.4. Results 

To infer the SFH of NGC 4395, we combined both fields 
to be able to directly compare with the cluster age distri- 
bution. The differences in our completeness limits for each 
band and each field are smaller than 0.1 magnitudes, al- 
lowing us to use the completeness values obtained in our 
first field (results for the second field can also be used in- 
stead) to study the SFH of the whole galaxy. In the SFH 
reconstruction, we used a resolution of 200 x 200 pixels for 
the Hess diagrams and a Gaussian kernel with a standard 
deviation of 0.02 mag in V — I (or B — V). Based on our 
observed Hess diagram we used 2 rectangular fitting areas 
{boxes2) for the SFH reconstruction with the same limits as 
described in the previous section. To check the consistency 
of our results we performed the analysis using both B — V 
and V — I colors, using the same parameters and boxes to 
fit the data. 
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Fig. 18. SFH for the galaxy NGC 4395. We studied the 
SFH for three different assumptions about the binaries try- 
ing to cover possible combinations between binary fraction 
and mass ratio for the colors V-I(left panel) and B-V(right 
panel) . 



Figure [T71 depicts the observed (left) and best-fit model 
(right) Hess diagrams (assuming / = 0.5 and q = [0.1, 0.9]). 
The final combination of parameters used to estimate 
the SFH were Padova 2008 isochrones, a Salpeter IMF 
(M=[0. 10,100] Mq), LMC metallicity, three different as- 
sumptions for the binarity, distance modulus of 28.1, a total 
extinction of Ab = 0.15, and the photometric errors and 
completeness limit obtained in Sect. 4. 

The phases that are most accurately reproduced by 
the fit are the main sequence, red core-He burning, and 
RGB/AGB phases. For the blue core- He burning phases, 
the model fit is less good, being redder than observed and 
with a more pronounced gap between the main sequence 
and blue core-He burning stars in the model than in the 
data. This is a generic problem with the isochrones that 
exists for any star formation history we adopt. 

Our most accurate estimates of the extinction were ob- 
tained by trial and error. We compared the main sequence 
stars from the fitted Hess diagram with data for those ob- 
served and assumed different values {Ab = [0.1,0.15,0.2]) 
for the total extinction until we obtained the best fit be- 
tween the two Hess diagrams (considering the main se- 
quence only). The difference between the foreground value 
{Ab = 0.074) and our best-fit model value {Ab ^ 0.15) sug- 
gests a low internal extinction in NGC 4395 of Ab ~ 0.076 
mag. This estimation was performed for LMC metallicity. 

The estimated SFHs (for V — I and B — V) are shown 
in Fig. [181 As discussed in the previous section, the SFHs 
become very uncertain at ages greater than a few hundred 
Myr; the apparent burst at ~ lO^ '' yr is probably not real. 
There is also a hint of a rapid drop at very young ages in V-I 
(~ 10^'^ years), but this may be caused by uncertainties in 
the isochrones. Furthermore, the youngest age included in 



Table 3. Average star formation rate [in MoYr^^] over the 
past 10^ years, estimated using different binarity, colors, 
and isochrones. 
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our artificial Hess diagrams is 4 Myr; if younger stars were 
present in the field, these would be included in the youngest 
bin. The appropriate lower age limit depends on how long 
young stars are embedded in their native molecular clouds. 
We present the SFHs for different binarity assumptions be- 
cause this parameter is also uncertain. We use the same 3 
choices in Fig. [18] as those used in Fig. [T3| above. Taking 
the values for each bin (in linear units), we estimated the 
average star formation rate for ages between 10^ and 10* 
years using the colors V-I and B-V, as shown in Table [3l 

There is an apparent increase in the SFR between 
logr - 10* ° and logr ~ 10'^ °, of ~ 0.5 dex. Based on 
Fig. [THl the average star formation rate over the past 10* 
years is SFR w 0.17 Mgyr"! within the two ACS fields 
(assuming / = 0.5 and q — [0.1,0.9]). We note that this is 
a factor of 4 higher than the global SFR for NGC 4395 that 
follows from using the far -infrared luminosity described 
in iLarsen fc Richtleij (|2000[ ). who used the calibration by 
iBuat fc Xul (|1996[ ). This highlights the difficulty in deter- 
mining extragalactic star formation rates from integrated 
light. 



5.5. Ages and masses of clusters 

After identifying the star cluster candidates (see Sect. 4) 
in the two fields, we need to estimate their masses and 
ages. These parameters w ere estimated using the pro- 
gram AnalySED, created bv I Anders et all (|2004l ). This pro- 
gram determines the ma sses and a ges using GALEV SSP 
evolutionary models (Schulz et al.l i2002'). The parameters 
used for AnalySED are a Salpeter IMF (jSalpeted [19551) 
i n the mass range 0.10 to 100 M0, Padova isochrones 
(jGirardi et al.ll2000[ ). and a LMC metalhcity. AnalySED 
performs a match, comparing the observed and theoreti- 
cal spectral energy distributions (SEDs), based on GALEV 
SSP models. As output, AnalySED provides, among other 
parameters, an estimation of the mass and the age for each 
cluster. 

In addition to our magnitude limit of = 23 {My ^ 
—5) for our clusters, we also define a lower mass limit 
at 10^'^ Mq to increase the likelihood that extended ob- 
jects are true clusters and not just chance projections of 
a few bright stars. However, over most of the age range 
{Log{age) > 7), these false detections would fall below our 
magnitude limit. 

Figure [TOj shows the age-mass diagram of the cluster 
candidates in NGC 4395. The number of clusters detected 
is rather small because the two HST pointings do not cover 
the whole galaxy and our coverage for clusters is reduced 
further by our requirement of the WFPC2 pointings for age- 
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Fig. 19. Age-mass distribution of the clusters detected in 
our two fields observed. The red dashed line represents 
My = — 5. Symbols are the same as in Fig. [51 



dating. We see that neither massive (M > 10^ M©) nor old 
(r > 10^ yr ) clusters are detected, in agreement with pre- 
vious work (|Larsen k Richtleijll999l: iMora et al 2009). We 
identified few clusters wi th t > 100 Myr for LMC metal- 
hcity in agreement with iMora et al.l ([2009) . although the 
number of objects was greater in their work. Nevertheless, 
the number of objects detected at these ages (and greater) 
varies, as can be see in Fig. 7 of Mora et al., depending on 
the metallicity used to estimate the parameters. 



6. Summary, discussion, and conclusions 

We have described our procedures to obtain the age dis- 
tributions of field stars and star clusters in HST images of 
nearby (13 « 4 Mpc) galaxies, in particular, we have de- 
scribed our implementation of the synthetic CMD method 
and tests of our program developed for this purpose. We 
found that the code recovers the star formation histories 
of synthetic populations with good accuracy, while errors 
in the derived SFHs of true populations are dominated by 
systematics errors. We have derived the SFH for NGC 4395 
for three different metallicities and find an approximately 
constant SFR over the past few hundred Myr assuming 
LMC-like metallicity. Only a modest amount of internal ex- 
tinction {Ab ~ 0.08) in NGC 4395 is required in addition 
to the Galactic foreground extinction to match the data. 
Hence, uncertainties in the total extinction is not a major 
source of error compared to, for example, binary stars and 
the choice of isochrones. Different assumptions may lead to 
changes in our estimated SFRs of up to a factor of 2-4 in 
specific age bins, although the global average will be less 
affected. Poissonian errors do not contribute substantially 
to the uncertainties because of the large number of field 
stars used in our analysis. 



After estimating the SFH of the field stars and both the 
age and mass distributions of star clusters, we can now mea- 
sure which fraction of stars are in clusters. We compare the 
fractions for cluster ages between 10^ and 10® years, since 
disruption due to tidal shocks (GMC and spiral arm en- 
counters) and evaporation is expected to become significant 
at older ages, while younger objects may still be unbound 
and prone to "infant mortality" . 

A rough estimate of the cluster formation rate (CFR) 
can be obtained by dividing the total amount of mass in 
clusters observed by the age range. From 10^ < t < 10® yr, 
the total mass observed in Fig. [TH]is Mobs ~ 8.3 x 10^ M©, 
which infers a log(CFi?) w —3.03 [M0yr~^]. We correct 
this for the smaller area covered by the WFPC2 camera 
compared to the ACS (a factor of ~ 2.27) and thus obtain 
log(CFi?) PS —2.68 [MQYr~^] for the same equivalent area 
covered by our field star data. 

A magnitude-limited sample does not allow us to ob- 
serve all the clusters in the age range. Using the observed 
total mass and assuming a cluster mass function (^(m)), 
we can es timate t he to tal mass in clusters for a certain 
age range. [LarsenI ( 20091 ) showed that a Schechter function 
(with = 2 X lO^M©) can model the initial cluster mass 
function in present-day spiral discs. However, the canonical 
cluster mass function, a power law with index —2, can be 
used for the same purpose. To estimate the cluster forma- 
tion rate we used both functions and calculated the total 
mass in the age range mentioned above. Since the most 
massive clusters observed have M< 10^ Mq, it makes little 
difference whether we adopt the Schechter function or an 
untruncated power-law. 

The total mass in the cluster system was estimated us- 
ing the equation 



{t2-ti) XTj m^{m)dm 
rt2 j-m^p rn,^ (m) dm, dt 

Jti ^nii^^u) ^ ' 



(2) 



where ii,t2 represent the time interval (t\ = 10^ years, 
h = 10® years), mio and TfXiip &-r6 the lower and upper mass 
limits over which the mass is normalized {mio = 10 and 
TUup = 10^ M0), and mnjy^i^i-^ is the (age- dependent) detec- 
tion limit of our cluster sample (dashed line in Fig. [TO]) . 
Since the last parameter, mum, is derived from the SSP 
models used, this integral must be evaluated numerically. 
Equation ^ assumes a uniform cluster age distribution, 
which is the most conservative choice we can make given 
the small number of clusters in our sample and the uncer- 
tainties in the age determinations (Fig. [B]). We also assume 
that 5'(r7i) is age-independent, meaning that we ignore dis- 
ruption. Applying the correction for undetected clusters, 
we obtained values for the CFR of log{CFR) = -2.37 
and log(Ci^i?) — —2.35 for power law and Schechter func- 
tions respectively, values that are corrected for the area 
covered by the HST detectors mentioned above. We have 
also tried applying Eq. ([2]) to 3 sub-bins in age, which 
changed the total overall CFR by a very small amount. 
We note that, strictly speaking, we are not measuring the 
true cluster formation rate, but rather the current age dis- 
tribution of surviving clusters. Dividing these numbers by 
the SFR {SFR = 0.17 MgYr-i), we obtain a ratio of 
F = CFR/ SFR = 0.0263. 

Our es timate of F is lower tha n the v alue F ^ 0.08 ob- 
tained by iBastianI (|2008D . iGielesI (|2009D estimated F for 
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Fig. 20. Relation between the cluster formation efficiency 
(r) and the total formation rate {TFR = SFR + CFR). 
An apparent power law function fits the date following the 
equati o n F » TFR^-^ . Symbols represent: ★ are data from 
iGielesI ([20091) for the galaxies M74 MlOl and M51; A is 
for the SMC (|Gieles fc Bastianll2008D : and o our estimation 
for NGC 4395. 



three galaxies namely M74, MlOl, and M51. The clus- 
ter formation rate was estimated by comparing theoretical 
with empirical luminosity functions. Giel es suggested that 
r has the tendency to increase with SFR. iGieles fc BastianI 
(|2008f) estimated F w 0.02 for the SMC using an accu- 
rate determination of the mean global star formation rate 
bas ed on the field st a rs (SF Rpa 0.1 MoYr"^) performed 
by [Harris fc Zaritskvl (|2004f ) and the CFR calculated in 
their work. Other st udies have found the same value for 
Goddard et aI1l2010t) . We took the 



F in the SMC (see 



(see Table 1 in his article), and 
and plotted the total formation 



data from iGieles (|200' 
IGieles fc BastianI (|200i 
rate {TFR = SFR+CFR) against F for these four galaxies 
and included our CFR and SFR for NGC 4395 to attempt 
to detect any correlation. Figure BOl shows the plot of log(F) 
versus log{TFR). There is an apparent correlation between 
these two values, which can be approximated by a power 
law F (X TFR^-^ . Although only five points are plotted and 
the TFR in NGC 4395 does not cover the whole galaxy, 
Fig. 1201 does suggest that the relatively low F value we de- 
rive for NGC 4395 is consistent with a general trend. In a 
future study (Silva- Villa et al. 2010, in prep.), we wiU in- 
clude more points in this plot for galaxies spanning a range 
in TFR and check the validity of this result. 

[Goddard et al.l (|2010l ) found an empirical relation be- 
tween F and the SFR density (EsFfi [Mgyr-^Kpc"^]). 
We estimated Y^sfr^ to be - 4.65 x 10"^ MeYr^^Kpc^^ 
for NGC 4395 whic h leads to F « 8% based on the 
I Goddard et al.l (|2010[ ) relation. This value is higher than 
that estimated in this work F ?» 3%. We point out that the 
area covered by our observations is not the total area of the 



galaxy and so the SFR density found in this work might be 
higher than the global average. 
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